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Abstract — In the target tracking and its engineering applica- 
tions, recursive state estimation of tlie target is of fundamental 
importance. Tliis paper presents a recursive performance bound 
for dynamic estimation and filtering problem, in the framework 
of the finite set statistics for the first time. The number of 
tracking algorithms with set-valued observations and state of 
targets is increased sharply recently. Nevertheless, the bound 
for these algorithms has not been fully discussed. Treating 
the measurement as set, this bound can be applied when the 
probability of detection is less than unity. Moreover, the state is 
treated as set, which is singleton or empty with certain probability 
and accounts for the appearance and the disappearance of the 
targets. When the existence of the target state is certain, our 
bound is as same as the most accurate results of the bound 
with probability of detection is less than unity in the framework 
of random vector statistics. When the uncertainty is taken into 
account, both linear and non-linear applications are presented to 
confirm the theory and reveal this bound is more general than 
previous bounds in the framework of random vector statistics.In 
fact, the collection of such measurements could be treated as a 
random finite set (RFS). 

Index Terms — lEEEtran, journal, BTgX, paper, template. 



I. Introduction 

IN the target tracking and its engineering applications, 
recursive state estimation of the target is of fundamental 
importance [IJ. However, the tracking system may not receive 
the information of the target, which is result from the probabil- 
ity of detection < 1 . Moreover, even though a measurement 
is received by the tracking system, it is hard to determine 
whether it is produced by the target or not, when the false 
alarm Pfa > 0. Therefore, at any time step, the number of 
measurement is random and it is unlikely to know whether 
there is missing detection or there are false measurements ||2l, 
0. 

This paper presents a recursive performance bound of 
estimation error with set-valued observations and state of 
targets, which is more general than the one with random vector 
measurement and state. The estimation problem where the both 
the measurement and the state are finite set, is very important 
in defense and surveillance [2|, [3|. The reason is that we 
cannot determine whether the target exists or not from the 
measurement and meanwhile the existence of the target varies 
with the time passing. In fact, the bound in the framework of 
random vector statistics is only a special case of our bound, 
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when the existence of the target state is certain. Therefore, 
this bound is a limit of a dynamic estimation error for the 
problem that the state set of the target is a Markov process, 
the measurement set is statistically dependent on the existence 
of the target, and the number of the points in the state set as 
well as measurement set is random at any time step. 

The error in this paper is a distance between the state set and 
the estimation set and thus the usual definition of Euclidean 
distance error for the random vectors cannot be applied. To 
solve this problem, a distance named Optimal Sub-pattern 
Assignment (OSPA) is given in H. OSPA is widely used in 
the performance analysis of algorithms (e.g. ||6] and 0), in 
the framework of the finite set statistics. 

Based on the OSPA, a mean square error (MSE) between the 
state set and estimation set is defined in this paper We want 
to find a limit of this MSE. When the state, measurement, 
and estimation are all random vectors, the limit for MSE 
is called Posterior Cramer-Rao bounds (PCRLB) In the 
framework of random vector statistics, is the measurement 
vector depends on the state vector x/j at time step k. The 
estimation vector icfe is based on the information gotten from 
the measurement before time step fc + 1: zi, Z2, • ■ ■ , z^. Cor- 
respondingly, in the framework of the finite set, the estimation 
set Xk is a function of all measurement sets: 



(1) 



Therefore, the definition of MSE between the state set and 
estimation set relates to the serial measurement sets as in ([T]). 

When the bound of the MSE is deduced, the PCRLB is 
also used. The PCRLB in |71 is a fundamental contribution 
for the development of the PCRLB. As the developments of 
the PCRLB in |7|, in the case of P<j < 1 and PpA = 0, 
Information Reduction Factor (IRF) PCRLB ||8| and enumer- 
ation (ENUM) PCRLB |9| can be applied by considering the 
effect of uncertainty in the measurement origin. Recently, these 
bounds are further tightened in [lO] and [11] respectively in 
clutter environment. Comparing to other PCRLBs, the ENUM 
PCRLB is the most accurate and the true bound for the case 
of Prf < 1 and - flU. 

There is another error bound based on OSPA given in (lT3l 
recently. This bound has a great influence on the derivation 
of the error bound in this paper However, the bound in llT3l 
models the state as a random set X, which is not a function 
of time step fc. In the other words, the bound in |13] is not 
recursive. Obviously, the meaning of a non-recursive bound is 
limited to the tracking system. 
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In this paper, for dynamic estimation and filtering problem, 
the state set Xk is a Markov process. In order to discuss the 
appearance and disappearance of the targets, the Xk may be 
{xfc} or empty according to the probability. At the time step 
k, the measurement set is: 



(2) 



Moreover, dU is modeled in the case of Pd < 1, whose 
influence is significant to the calculation of error bounds. 

In addition, part of this result in this paper has been reported 
in ifTSl . where only linear filtering case was presented and the 
discussion of the results was absent. 

Section revises the traditional PCRLB and the basic 
knowledge of random set statistics. A new concept of mean 
square error (MSE) cr^ between Xk and Xk is defined in the 
section Hill In the section IIVI a recursive form of this error 
bound is derived. This bound is discussed in section |V] We 
present two numerical examples in the section |VI] Proofs of 
the propositions are in the Section IVlIl Conclusions are drawn 
in the Section IVIIII 

II. Background 

A. Recursive Form of the PCRLB 

For a random vectors filtering problem, the state dynamic 
equation is given by: 



Xfe+i = ik (xfe, Wfc) 



(3) 



where is the state transition function, and is a zero- 
mean white Gaussian process noise, with covariance matrix 

When the target is detected, the measurement equation is 
given by: 

Zfe = hfc (xfe, Vfe) . (4) 

where is the observation function, and is a zero-mean 
white Gaussian noise, with covariance matrix R/j. 

Let Xfe be an unbiased state estimator based on the sequence 
of measurements {zi, • • • , z^}. The covariance of this estima- 
tor has a lower bound expressed as follows |[l]: 



E 



(Xfe - Xfe) (Xfe - Xfc) 



(5) 



where is referred to as the Fisher information matrix (FIM), 
and the = J^:^ is the PCRLB. 

As in f9^|, when the target is detected, the recursive formula 
of FIM is as follow: 



'fc+i 



Also as in f9l, if the target is missed, or the FIM is for 
predictive, the recursive formula of FIM reads 



Hfe 



Vxjhfe (Xfc)]' 



(7) 



Jfe+i - Qfc ^ 

- Q-^E {Yk} [Jfe + E {FjQ-iFfe}] E {F,.} Q-^ 

(8) 

In a word, whether or not z^ ever exists, the FIM at fc+1 can 
be calculated by the dynamic equation, measurement equation 
and the FIM at k. 

B. Random Finite Set 

Random finite set (RFS) is a random variable which takes 
value as finite set lfT3l . The element of this set is unordered 
random variable and the number of the elements is random 
and finite. Finite set statistics (FISST) is developed by Mahler 
||3| and widely considered an effective tool for the multi- 
target tracking system. In the perspective of modeling the 
tracking system, two types of RFS are often used: Poisson 
and Bernoulli RFS. Based on the model of Poisson RFS, a 
filter named Probability Hypothesis Density (PHD) filter [T?] 
are applied in several fields [16] , I.17J . However, PHD is a first- 
order statistical moment of the multi-target posterior [lT4l . and 
Poisson RFS is apt to model the multi-target tracking system. 
Therefore, Poisson RFS model does not suit to the single target 
appearance and disappearance problem in this paper The filter 
derived from Bernoulli RFS attracts substantial interest and is 
used widely recently [14J, 1 18,| . As in [14], here a Bernoulli 
RFS on a space S is defined by two parameters r and p (•): 




X = {x} ; 
otherwise. 



(9) 



where the f {X) is the density of the RFS X on the space of 
finite sets. 

For the function g taking value on the set X, the set integral 
of this function is IS): 

/ g{X)5X = g{fb)+Y.-^ I 5 ({xi, • • • , x„})dxi, • • • , dx„ 



(10) 

The expectation of the function g on a RFS of density / is 



E[h] 



h{X)f{X)SX 



(11) 



If the state of the target is X, the estimation is X [Z), 
where Z is the measurement of the target. The distance defined 
between two sets X and X [Z) is as follow |13J: 



Qll'E {F,} [J, + E {F^Q^T^F,}] ^E {F^} ^-^^ 

(6) 

where the matrices F^ and are respectively the Jacobians 
of nonlinear functions ik and h^: 



e(^X = {x},X(Z) = {x}j=x-x; (12) 
e = 0, 1 {Z) = {x'}) = Go, for any x'; (13) 
e(x = {x},X(Z) = 0) =ei,/or any x; (14) 

e(x = 0,l(Z) = 0) =0. (15) 

Since the number of element in the set may be zero or 
one, the difference between X and X (Z) is defined in (flj] ) 
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and ( fT4l i. when the numbers of element of this two sets are 
different. For one thing, if there is no target in reality, we still 
estimate there is one target, the error is Gq. Or perhaps, there 
is one target, but we estimate there is no one, such error is ei. 
In a word, Gq and ei indicate the mismatches of cardinality. 

According to the definition of the error, the mean square 
error between X and X [Z) is given as in ifTSl : 



11 = E 



e[x,X{Z))e(x,X{Z)y 



= / ^ e(x,X{Z))eix,X{Z)) p{X,Z)SXSZ, 



(16) 



where p {X, Z) is the joint density of the state set X and 
measurement set Z. 

III. Necessary Definition 
A. Observation-sets Sequence 

In the framework of random vector statistics, at time step 
fc, the estimation is a function of all measurements from time 
step 1 to fc, and contain all information of such measurements. 
These measurements are Zi, Z2, • • • , z^, and thus the estima- 
tion is Xfe ({zi, • • • ,Zfe}). Correspondingly, in the framework 
of random set, firstly, we should determine how to state the 
measurement sets from time step 1 to fc, which is defined as 
observation-sets sequence Qk,n- 

At time-step fc, the possible time-sequence of observation- 
sets is given as follows: 



(17) 



where „ denotes the measurement is empty or not for 
sequence number n, at time-step fc, and thus, n = 1,2, ■ ■ ■ ,2'^ . 

In order to not only simplify the form of the bound deduced 
in this paper, but also indicate the influence of whether the 
observation is empty or not, the arrangement of the elements 
in &k,n is not random, but follows the rule: 

When k = l 



5,ei,2 = {zi} 



When fc = 2 

62,1 ={0,0}, ©2.2 = {{Z1},0}, 

62,3 = {0,{Z2}}, ©2,4 = {{Z1},{Z2}} 

When k>2 

0fe+l,ri 

{©fc,„,0}, l<n<2'= 
{6fc.„_2'^,{zfe}}, 2'^- + 1 < n < 2^-+i 

{efc_i,„,0,0}, 
{6fc_i,„_2fc-i,{2fe},0} , 

{0fc-l,n-2'',0, {^fe}} , 
6fc-l,„-2'=-2'»-i 
{Zk},{zk+l} 



(18) 



(19) 



1 < n < 2'=-! 
2k-i + i<n<2'' 
2*= + 1 < n < 2*= + 2'=-! 



set is empty Z^+i^n = 0, where the sequence number n is in 
the scope 1 < n < 2*^'. This situation appears when there is 
no target or the target is missed. The second part is in the 
condition that Zk+i.n 7^ 0, where n is in the range 2^^ + 1 < 
n < 2^^+^. This condition results from that there is a target 
and it has been observed. 

It is notable that, at time step fc + 1, the 6^ „_2'» with 2^ + 
1 < n < 2^^+^ and the Qk.n with 1 < n < 2''' share the same 
observation-sets sequence, which is all possible observation- 
sets sequence at time step fc. 

B. Error Bounds Defined based on Observation-sets Sequence 

For certain time-sequence of observation-sets 6fc „, the 
estimation can be written into a particular form. Then, the 
estimated error is defined as following: 

= /•■•/ CkMXk,Q,^n)6Xk6Z,,^ ■ . . <5Z,,„ > P,,„ 

(21) 

where 

Ck,n = e (Xfc,Xfc(6fc,„)) e(Xfe,Xfe(6fc,„))^ (22) 

At time step fc, Pk,n is the error bound for the particular 
observation-sets Qk,n- 

Because Qk,n can cover all possible conditions of observa- 
tion when the number of sequence n takes value from 1 to 2^, 
the total error between Xk and its estimation Xk{Zi ■ ■ ■ Zk) 
is Sfc,„: 



J... J e(^Xk,Xk{Zi---Zk))e[^Xk,Xk{Zi---Zk)) * 
p{Xk,Zi, ■■■ , Zk)SXkSZi ■ ■ ■ SZk 

= £/ •■•/ ''{^k,Xk{Qk,n))e(^Xk,Xk{^k^n)Y* 

n=l p{Xk, 6fc,„)5Xfc5Zi,„ • • • SZk,n 

2'° 

(23) 

Therefore, the error bound of such total error is as follow: 



Efc = ^ I]fc,„ > Pfc - ^ Pfc, 



(24) 



ri = l 



Hence, in the rest of paper, what we deduce is the recursive 
form of Pk,n- According to (l24l l. the total bound Pfe can be 
obtained by the sum of 'Pk,n- 

IV. Recursive Form of the Bound 
A. Random Finite Set Models 



Since the target in either 'present' or 'absent' state, the state 
2fc-i + 1 < 7T, < 2^^+^ of the target is modeled by Bernoulli RFS mentioned in section 



(20) 

(I20I 1 shows that the elements in &k+i,n can be divided into 
four parts with equal number of elements, according to the 
measurements at time step fc and fc + 1. If the dividing is just 
on the basis of the measurement at time step fc + 1, there are 
two parts. The first one is in the situation that the measurement 



El 

For the dynamical model, the Markov transition density is 
defined by: 



f{Xk+i\Xk = {^k}) = 



rV'(xfe+i |xfc), Xfe+i = {xA;+i} 



l-r. 



Xk+i — 



(25) 
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r(Y IV a\ / (l-OPo(Xfe+l), Xfe+i={Xfc+i} 

(26) 

where r G [0,1] represents the probabiHty of the state of the 
target at time step fc + 1 survives from the state at time-step k 
or remains empty. It means, conditional upon Xk = {xfc}, this 
target disappear with the probabiHty 1 — r. If there is no target 
at time-step , a new target would bear with the probability 1— r, 
and the initial density is pq (x^+i). Therefore, r is named as 
the maintenance probability. ■0(xfc+i ) is the probability 
density of a transition from xj. state to x^+i. 

The prior probability function of the state set is also 
Bernoulli RFS: 



bpo (xo) ,Xo = {xo} 
1 - &, = 



(27) 



where b e [0,1] representing the probability of the target 
existing initially, b is named as the initial probability. 

The probability of detection is < 1. The measurement 
model is: 



g{Zk\Xk^{^k}) 



g{Zk\Xk^(/}) 



Pdii'Z'k \^k),Zk = {zfe} 
l-Pd,Zfe = 

0,Zfc = {zfe} 

i,Zk = 9 



(28) 



(29) 



where ^(zj.|xfc) is the measurement likelihood when the 
target is existing and detected. ( |28] | indicates there is some 
uncertainty in detection, and (|29] l means there is no false 
observation. 

B. Derivation of the Bound 

1) Derivation ofPk- As discussion in section HiH in order 
to calculate the recursive form of Pfe, we should get the 
recursive Pk,n- 

Proposition I: When time-step fc > 0, n = 1, • • • , 2*^+^ , the 
bounds sequence Pk+i,n obeys the recursion: 



Pfe+l^n' 1 < n < 2^= 

Pfc+1,,1 = { [Jfc+i,„] ^ * Pr (efc,„_2fc , Zk+1 ^ 0) , 

2*^ + 1 < n < 2'^'+! 



otherwise 



(30) 



(31) 



where 



n+i,n = e^ef (Pr (Gfc,™, Zk+i = 0) - Pk+i,n) (32) 



Pfe+i,™ = eoeJpfe+i,„ + [Jfc+i,n] ^ * Pr (6fc,n, Zk+i = 



(33) 



Pk+l 
A 



I ■ ■ ■ I Pi^k+l = 0, 0fc,n, Zk+l,n = 9)6Zi.n ■ ■ ■ SZk,n 

(34) 



Pk+i.n indicates the probability of the state set of the target 
is empty at time step fc + 1, when all the measurements from 
time step 1 to fe + 1 are all known. 

The proof of this proposition is in the section I VII I 
From proposition I , we can see that the problem of 
recursion of Pfc.„ reduces to the recursion of Jk+i,n, Pk+i.n 
and Pr (Ofc+i.n)- Then we deduce how to obtain the recurrent 
formulas of all these three factors. First, based on (|6) and dH), 
for a particular Qk.n defined in (|20] i. the FIM ik+i,n can be 
calculated from „ when fc > 0: 



Jfc+l,n — 

' Qfc ' - ^k ^E {Pk} [Jfe.n + E {F^Q^T^F,.}] {Pk} Q 

1 < n < 2^= 

-Ql'E{Pk] [ik,n-2-+E{PlCl-^Pk)Y^E{Pk}Cll 

2^ + l<n< 2'=+! 

(35) 

If there is no process noise w^, then ( |35] | reduces to: 

[ [p-Y]':ik.nPl\l<n<2^ 
I Hj+iR^+iH,+i,2^- + l<n<2'^-+i 

(36) 

The initial FIM is calculated from the prior probability 
function po (^o): 



Jo,i -£^{-A^° logpo (xo)} 



(37) 



2) Derivation of Pr(9fe+i.„).- Secondly, Pr(9fe+i.„) 
consists two parts: Pr (9fe.„, ^fe+i = 0) and 
Pr (9fc „_2'«i ^fc+i 7^ 0), according to the Qk+i,n- When 
fc = 0, the initial probability of the empty measurement is as 
follow: 

Pr (9o,i, - 0) = Pr (9i,i) = 1 - & * Prf (38) 
Pr (9o,i, ^1 7^ 0) = Pr (9i,2) = b*Pd (39) 

When fc > 1, the recursion of Pr(9fc+i,„) is given in 
Proposition II: 

Proposition II: When fc > 1 

Pr(9fc+i,„) 

Pr(9fe,„)p(Zfe+i,„ -0|9fc,„), 

1 < n < 2'=-! 

Pr (9fc,„) (1 - rPa) , 

2k-i + 1 < „ < 2'^- 

Pr(Bfe.„_20 [^-p{Zk+l.n = |©fe,„-2'=)] ' 

2*^ + 1 < n < 2'^' + 2'^-! 

Pv{ek,n-2^) rPd, 

2'^' + 2'^-! + 1 < n < 2*^+1 

(40) 

The proposition II is proved in the section IVIII 
At time step fc + 1, for the number of queue n in the two 
range that 2''-'^ + 1 < n < 2^= and 2'' + 2''-^ + 1 <n< 2''+^, 
the measurement Zk.n 7^ 0- It means that, the target exists 
at time step fc, because it is assumed that there is no clutter 
Then, at time step k-nl, the state of the target can be written by 
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the state transition model. Hence, p {Zk+i.n 7^ |0fe,ri ) can be 
calculated easily. However, when the measurement Zk.n — 0, 
which corresponds the two range that 1 < n < 2*^^^ and 
2'° + 1 < n < 2*"' + 2*^^^, it is uncertain the reason is the state 
set is empty or there is a miss detection. Therefore, it is hard 
to determine p{Zk+i.n = \'S>k,n) or p(Zfc+i,„ ^ \Qk,n), 
when Zk^n = 0. 

The recursion of p {Zk+i.,i = 9 \Qk,n) is given in the 
proposition IV, for the recursive form of pk+i,n is pre- 
sented in the proposition III, which is also related to 

p{Zk+l,„ = ISfc^n)- 

3) Derivation of pk+i,n-' Pk+i,n is defined in (l34l i. When 
fc = 0, the initial pi^i is calculated as: 



= 1-6 



(41) 



When fc > 1, pk+i,n is calculated as the proposition III. 
Proposition III: When k > 1 



Pk+l,n — 

1 <n<2'' 

(42) 

The proposition III is testified in the section IVIII 
When the measurements from time step 1 to fc + 1 have 
all obtained, pk+i,n shows the probability of the state set is 
empty at time step k + 1,. Hence, it is obvious that pk+i,n is 
a function of p {Zk+i,n = \Qkjn )■ 

4) Derivation of p [Zk+i.n — |0fc.n)-' From the proposi- 
tion II and proposition III, we can see that both Pr (Ofe+i.„) 
and pk+i.n are functions of p(Zfc+i.„ — \Qk,n) when 1 < 
n < 2^ . Therefore, the key is how to get the recursion of 

p{Zk+X,n = |6fc,n)- 

Proposition TV: When fc > 1 



p(Zfe+i.„ = 0|efc,„) 

1 - rPd + Pd (2r - 1) * 

p(Zfc.„=0|efc_i,„)-(i-Prf) 



PdP(2fc,„=0|efc- 
1 - rPrf, 



1 < n < 2*^-1 



= r(p(Zfe,„ = 0|efc_i,„)) 



(43) 

The proposition IV is testified in the section IVIII 

By taking ( l43T l into the Proposition II and III, Pr (9fc+i^„) 

and pk+i.n can be obtained from which at the time step fc. 

Therefore, the Pfc+i.„ can be calculated in a recurrent form. 

From the (l24l i. the total error bound can be obtained. 



Combining the definition of pk+i.n in (O, setting Gi = Gq, 
it is easy to derive that 



/ • • • / /(Xfe+i = 0, efc,„, Zk+l,n = 0)<5Zi,„ • • . SZk.n > 
e^ej -[Jk + i f r ^ 0, 0fc,n, ^fc+l.n = 0) 

e,ef + [J. + i,„]-^ J ■■■J SZi^n ■ ■ ■ 5Zk,n 

(46) 

While the number of the scans of measurements increases, 
it should meet the condition that 



e^el = e^ef > [J/c+i.n] ^ 



(47) 



As a result: 



/ • • • / fi^k+l = 0, 0fc,ri, Zk+l.n = 9)5Zi^n ' ' ' SZk,n > 
J ■■ ■ J f{^k+l 7^ 0, 0fc,ri, Zk+l,n = ^)6Zi^n ' ' ' SZk,n 

(48) 

It denotes that, when the probability of empty state set is 
more than which of not empty, the bound is in the form that 

Pfc+i,n = ^1^1 (PrCQfe,"; = 0) - Pk+i.7i)- In the other 
word, there is Xk+iiQk,^ Zk+i,n — 9) — 9, when the bound 
is attained. 

On the other side, if ( |47] i is satisfied, and if 

/ • ■ • / fiXk+l / 0, efe,„, Zk+l,n - 0)<5^1,„ • • ■ SZk,n > 
I ■■ ■ I fi^k+l = 0, 0fc,n, •Z'fc+l.n = 9)SZi^n ' ' -SZk^n, 

(49) 

the estimation is Xk+ii'dk.n, Zk+i,n = 0) ^ 0, and 



the bound is that P 



k+l,n 



Gq^O Pk+l,n + [Jfe+1,; 



Pr (0fc.n, Zk+1 = 0). Then, this bound can be compared with 
the PCRLB. 



B. Comparison with Previous Results 

The enumeration PCRLB has been verified as the exact 
bound in the case of Pd < 1, both in a linear and a 
nonlinear case. It is the optimal lower bound for tracking in the 
framework of the finite vector statistics. Rewrite the PCRLB 
computed via enumeration in jjO) and lfT2l as following: 



V. Discussion 

A. The Meaning of This Bound 

If the state of target is one-dimensional, the bounds 
calculated with empty measurements dSTT l turn to: 



fc+l,n 



When P 



fe+l,ri 



^k+l,m^k+l,n — ^ k+l,n 



-ffc+1 n' i'^ means that 



eQelpk+i,n + [Jfc+i,n] ^ * Pr {Qk,n, Zk+l = 0) > 
e^ej (Pr (Ofc^n, Zk+i - 0) - Pk+i,n) 



(44) 



(45) 



Pk+i{ENUM) - 

[3k,n + E{¥lq,-'¥k]] '£;{Ffc}Q-^ 
xPr(efe,„,Zfe+i =0) 

Qfc ^ + ^ {^k+i^kli^k+i) 

-Q^ii?{Ffe}x 

[J,.„_2. +ii;{FjQ,TiFfc}] 'i?{Ffc}Q^i . 

xPr(efc.„_2'.-,^fc+i^0) 

(50) 

At the same time, if Wi\ and ( |49] l are satisfied, our bound 
is given by: 



n=2''+l 
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fc+1 



2" 
n=l 



RFS = 
eoeoPfe+i,n+ 

Qk' -Qk'E{Fk}x 
[jfe,„ + i?{F^Q-iF,}] ' 

[ xPr(efe,„,Zfc+i =0) 

Qfc ^ + {Hfc+iR-fc+iHfc+i} 



> + 



2fc + l 
n=2*' + l 



Pr(efc,„_2'=,^fe+i 7^0) 



(51) 

It is obvious that ( fSOl l is equal to ( |5T] i. in the condition that: 

= !■■■! fiXk+i = 0, efc,„, Zk+i,n = 0)'5^i,„ • • • sZk,n 
= 

(52) 

It means that Pfc+i(RFS) is equal to Pk+i{ENUM), when 
the target exists form the beginning to the end. 

VI. Examples 

In this section, the previous theory is illustrated by two study 
cases: the one related to a linear filtering model and the other 
one referring to a non-linear bearings-only tracking. 

A. Linear Filtering Case 

This section illustrates the application of previous theoret- 
ical result by a linear case with Gaussian noise. When the 
target exists: 

1 r~i 



Xk Xk Uk Uk 



Xfe = {x,}={[ 

The target is detected: 

Zk = {Zfc} = |zfc = [ Zy^k ]^| 

The target motion is modeled as a linear equation: 

Xfc+i = FfcXfc + Wfc, 
1 T 
10 
1 T 
1 
Wfc ^ N (0, Qfc) , 
T3/3 TV2 
TOO 
r3/3 
T 



(53) 



(54) 



(55) 



Qk 



In this case, the time interval T is 5s. The intensity of the 
process noise q = 10^^. 

The measurement equation is as follow: 



Zfc = HfcXfc+Vfe, Hfc 
Vfc -7V(0,Rfc), 



10 
10 
Rfc = diag (cr^,cr^) 



(56) 











1 ^0 ^ 









The variance of measurement noise is = a"}, ~ 25^. 



The initial target state standard variance is c,. = 100, c„ = 5, 
d the initial FIM is J^} = diag{c'^,cl,cf,cl). 
The errors in cardinality mismatches are: 



(57) 



Fig.l, 2 and 3 show the root-mean-square error (RMSE) 
bounds between two sets and in (a) x-position and (b) x- 
velocity for ten scans. Similar results can be also obtained 
for the y-axis of position and velocity. 

1) The Influence of r: In Fig.l, the unchanged parameters 
used in the computation of the bounds are Pd — 0.8 and b — 
1. Here, is the probability of the target existing initially. We 
compare the different bounds at various value of the parameter 
. r is the maintenance probability. In addition, there is the for 
case using ENUM method as in (for short RMSE (ENUM)) 
Q. When the target exist from the first time-step to the last one 
(& = l,r = 1 for the models of RMSE (RFS)), this situation 
is the same as in which the RMSE (ENUM) is calculated. 
Moreover, RMSE (ENUM) is the tme bound for the case of 
Pd < 1 and Pfa = 0. 

In Fig.l, when r — 1, the RMSE (RFS) is approximate 
to the RMSE (ENUM), and as the number of the scans of 
measurements increases, they become the same. The reason 
is that, as discussed in section |V] when it is not satisfied 
the condition that ej^e^ = e^e]^ > [Jfc,„]"\ RMSE (RFS) 
and RMSE (ENUM) could not be the same but just close. 
With the number of the scans of measurements increasing, it 
become satisfied that ej^e^ — GqC^ ;» [Jfe,n]~^, and then 
RMSE (RFS) and RMSE (ENUM) are similar This problem 
would be further discussed in the following passage. 

When ^ r < 1, by taking account to the uncertain of 
the existence of the target, the bounds enhances. In reality, it 
means that, if the target disappeared with the probability 0.1 
(r = 0.9, h = 1), the performance of estimation would not 
reach the RMSE (ENUM), but the RMSE (RFS). In the other 
word, the calculation of RMSE (ENUM) is overly optimistic, 
when there is uncertainty about the target existing or not. 

In addition, we confine ourselves to the matter under the 
situation ^ r < 1 in the following discussion, when the 
initial probability of existing target h — 1. Because the case 
of0<r<Cl,6=l means that the target being preset or not 
is changed dramatically with time step. This is unrealistic and 
with little significance to discuss for tracking system. 

2) The Influence of Pd: As shown in Fig. 2, The RMSE 
(RFS) and RMSE (ENUM) are compared in the case of Pd = 
0.7 and Pd = 0.9. The unchanged parameters are r = 0.9 and 
6 = 1 for the models of RMSE (RFS). This means that the 
target enters at the first step, and then disappeared with the 
probability of 0.1. 

Fig. 2 shows that the RMSE (RFS) is always larger than 
RMSE (ENUM), which is also illustrated in Fig. 1 . The reason 
is, when r < 1, the uncertainty of the existence of the 
target improves the bound of estimation. In Fig. 2, when the 
probability of detection is reduced from 0.9 to 0.7, the RMSEs 
calculated by both the two methods are increased. However, 
the influence of the target appearance or disappearance is more 
significant than the influence of the miss detection. Because 
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♦ RM.SE(RFS),r=I 




(a) Position in X-axis 



♦ RM.SE(RFS),r=l 




(b) Velocity in X-axis 
Fig. I. Comparisons of RMSEs witli different r (linear filtering) 





HIH RMSEI RFS),r=0.9,P^=0.9 
^RMSE(ENUM).P^=0.9 
-X - RMSE( RrS).r=0.9,P^=0.7 
RMSE(ENUM),P^=0.7 
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Number of Scans.k 



(a) Position in X-axis 





















H|HRMSE(RFS),r=0.9,P^=0.9 










^RMSE(ENUM).P^=0.9 
-X . RMSE( RFS ),r=0.9.P^=0.7 










RMSE(ENUM).P^=0.7 



















































(b) Velocity in X-axis 

Fig. 2. Comparisons of RMSEs with different probability of detection (Linear 
filtering) 

the error in cardinality mismatches has great effect in the 
calculation of the RMSE (RFS), which is paid no attention in 
the calculation of the RMSE (ENUM). Therefore, the RMSE 
(RFS) would be the true bound in the case of < 1 and 
Pfa ~ 0, if the targets disappeared with certain probability. 

3) The Influence of Cardinality Mismatches: In order to 
illustrate the discussion in|V] we reset the errors of mismatches 
of cardinality Gq and ei: 











ei = 




7 ^0 ^ 


2C^; 



(58) 




(a) Position in X-axis 




NuinbiT of Scaiis.k 



(b) Velocity in X-axis 

Fig. 3. Comparisons of RMSEs with different r with high error in cardinality 
mismatches (linear filtering) 



The other settings are as similar as what in the Fig.l. 
Comparing the Fig.l and the Fig. 3, when it is satisfied that 
Gj^ef" ~ GgGj^" 3> [3k,n\~'^, it is evident that the RMSE 
(RFS) and RMSE (ENUM) are always the same. As the 
enumeration PCRLB is verified as the exact bound in the 
case of Pd < 1, our method calculating tracking bound can 
always attain the optimal method, when errors g^ and Gq 
are set as (ISSl l. However, as in Fig. 3, overestimated error in 
cardinality mismatches leads to an unreasonable higher RMSE 
(RFS). Moreover, the overrating RMSE (RFS) reduce slowly 
with the number of the scans decreasing, which would be 
meaningless for tracking system to some extent. Therefore, 
the setting of the errors in cardinality mismatches in (l57i is 
more reasonable. As in figure. 1, the RMSE (RFS) and RMSE 
(ENUM) approach each other as the number of the scans 
of measurements increases, and then they become the same 
at last. In reality, it is usual that the errors bought by the 
mismatches in the number of targets between the true state 
and estimation, i.e. Gq and e^, are designed according to initial 
FIM. 



B. Nonlinear Filtering Case 

This example is as similar as the bearings-only tracking case 
in fl2\. This system can be applied in electro-magnetic (EM) 
equipment, electronic warfare devices (ESM) and passive 
sonar llT2l .subsection text here. 

The observer, named ownship, is a moving platform car- 
rying sensor Its state vector is denoted as and assumed 
known. The target vector is denoted as x|,. The relative state 
vector is defined as: 

Xfc = - = [ Xfe Xfe 7fe 7fe ] ^ (59) 
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-ownship 
►-target 




Fig. 4. Bearing-only tracking scenario 



where (xfej7fe) is the relative target position and (xfe,7fc) is 
its velocity. The dynamic equation is as following: 



Xfc+i — FkXk — Uk,k+i 



(60) 



where is denoted in dSSl l, and the effect of a mismatch 
between the observer and the target motion model is accounted 
by: 



U, 



xl+^-Xl-Txl 
Xk+i ~ Xk 

il+i-il-TYk 
Yk+i - il 

The measurement equation is: 



Zfc = hk (xfe) + Vk 



where 



hk (xfc) — arctan 



Xk_ 

Ik 



(61) 



(62) 



(63) 



Vk is a zero-mean white with co variance Rk = cr^ = (1°) 



Based on (|7|i, the Jacobian of hk (xfc) is calculated as: 



Hk = 











(64) 



The initial target state standard variance is c,. = 
10000to,Ci, = lOOm/s, and the initial FIM is 3^1 = 
diag{cl,cl,cl,cl). 

Ownship is moving as a uniform circular motion. The 
angular velocity is w = 1.0125°/s. The dynamic equation 
of the observer is given by: 



(1 



sin(a;T)/a; 

cos(wT) 
- cos(wT))/w 

sin(wT) 



= e 



(65) 



(-1 + 



cos((a;T))/tj 
— sin(a;T) 

cos(cjT) 



(66) 



initial target 
150m/s 20fcm 



The 

[ -25fcm 
the initial observer 
[ -30fcm 200m/s 50fcm 



vector 

1 T 



vector 
r. The 



and 



target and 



state 

lOOm/s 
state 

Om/s 

observer trajectories are shown in Fig.4. 

Fig. 5 and Fig. 6 show the RMSE bound between two sets 
Xk and Xk{Zi ■ ■ ■ Zk) in (a) y-position and (b) y-velocity for 
twenty scans. Similar results can be also obtained for the x- 
axis of position and velocity. 



1) The Influence of r: In Fig. 5, the unchanged parameters 
used in the computation of the bounds are and . We concentrate 
on the influence of various value of the parameter . Since there 
is no process noise for both ownship and target, the calculation 
of FIM is based on ( l36b . 

As shown in Fig. 5, b — l,r — 1 for the models of RMSE 
(RFS) mean that the probability of target existence is unity. 
This condition can be compare with that of calculating the 
RMSE (ENUM) for they are the same scenario. 

Note that, in this bearings-only tracking case (Fig. 5), the 
RMSE decrease more steeply than which in the linear case 
(Fig.l). The reason is, at the initial several scans, initial FIM 
J(7j^ impacts the RMSE most, because the measurement of 
target is missed initially with high probability, comparing 
to the following scans. However, the covariance matrix of 
measurement noise Rfc influences the estimated error more 
with more measurements observed. In the Fig.5, there is the 
relationship that Jq^|^ R^, while in the Fig.l, the Jq^|^ is 
bigger than Rfe, but closely. Therefore, the bounds shown in 
Fig.5 is influenced by the initial FIM Jq \ more than which 
in Fig.l. 

Moreover, in the Fig.5, when r < 1, the bounds intersect 
the bound of r = 1. While, in the Fig.l, these bounds always 
exceed the bound of r = 1. The reasons are not only the 
initial FIM Jq \ ^ Rfe discussed above, but also the error 

e(^X ^ (1},X{Z) = 0^ =0, which is defined in O- The 
initial target state standard variance is Cr = 100m in the 
linear filtering case, while which in the nonlinear one is 
Cr — 10000m. Therefore, the estimation of absence of target, 
where e (^X ~ 9,X (Z) = 0^ =0, plays more roles in the 
bounds when r < 1 in the Fig.5, and contributes that the 
bound of r < 1 is lower than that of r = 1. In fact, by the 
accumulation of several scans, the bounds r < 1 must be no 
lower than that of r = 1, which is also shown in the Fig.5, 
because it take the uncertainty of the existence of target. 

2 ) The Influence of b: In order to indicate the influence of 
the target existing or not initially, the value of b is changed, b 
is the probabiUty of the target existing initially. In all previous 
examples, we consider the situations that target exist at the first 
step (b = 1). In the Fig. 6, we reset b = 0.1, which means that 
the target appear on the probability of 0.1 at the first time step. 
In this case, the bounds of RMSE (RFS) is hard to compare 
with the RMSE (ENUM), because when the RMSE (ENUM) 
is calculated, it is the same situation that b ~ l,r = Ifor 
the calculation of RMSE (RFS). Nevertheless, the relationship 
between RMSE (ENUM) and RMSE (RFS) is discussed in 
detail in the linear example and the need of comparing them 
is little in this example. Here we contain b ~ 0.1 and vary r in 
the Fig. 6. In Fig. 6, the probability of detection is set P^j = 0.9 
for all bounds. 

As shown in Fig. 6, the influence of 
e = 0,X (Z) = 0^ = is more significant than 
Fig.5. Because this case means that at first step, the state 
set of the target is empty with the probability 0.9, and 
then the cardinality turn to one with the probability 0.1 
(r — 0.9, b — 0.1). In the other word, the target enter with 
the probably 1 - r. When r < 1, the bound of RMSE (RFS) 
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(b) Velocity in Y-axis 
Fig. 5. Comparisons of RMSEs with diiferent r (nonlinear filtering) 
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(b) Velocity in Y-axis 
Fig. 6. Comparisons of RMSEs with b=0. 1 and varying r (nonlinear filtering) 



is bigger than RMSE (ENUM) in Fig.6, except the initial 
several scans. The optimality of the RMSE (RES) is verified 
again. It is noted that, in the case of , the bound of RMSE 
(RES) with r = 1 is always less than RMSE (ENUM). This 
means that the target births at the first time step with the 
probability 0.1, and keep the similar state as the first scan. As 
in Eig.6, to the situation that the target is rare, the estimation 
of empty set can reduce the error 



VII. Mathematic Proofs 

A. Proof of Proposition 1 

When Zk+i ^ 0, the error bound at time-step /c + 1 is as 
following: 

SXk+lSZi^n - ■ -SZk+l^n 
- Xfc+i) (Xfe+i - Xfe+i)^* 
= /•••/ P{Xk+l,n = {Xfc+l} , 9fe,„_2'=7 ^fe+l,ri = {Zfe+l}) 

dxk+iSZi^n ■ ■ ■ SZk,,idzk+i 
Qfc +^{Hfc+iRfe+iHj.^J - 
Q^'E{Fk}[Jk + E{FlQ^'Fk}r'E{F,}Q^' 

*Pr(efe+i.„) 

= [Jfc+l,n]"'*Pr(e,,„_2^,Zfc+i ^0) 

(67) 

When Zk+i = 0, the error bound is given follow: 



> 



= /•■•/ 



/•••/ 



Ck+l,n * f{Xk+l, 'dk,n, 
6Xk+lSZi,n ■ ■ ■ SZk.n 

e (^9,Xk+l{Qk,n, Zk+l.n = 

e(0,lfe+i(efc,„,Zfe+i,„ ==0)) * 

p{Xk+l = 0, Oa_.,„, Zk+l.n = 9)SZi^n ■ ■ ■ SZk,n 
e (^Xk+1 = {Xfc+i} , Xk+l{Ok,n, Zk+l.n = 0)) * 

+ /■■■/ ^i^^k+1^ {y^k+l} ,Xk+l{Qk,n, Zk+l.n ^'^)^ * 
/ {Xk+1 — {x/j+i} , 6fc,„, Zk+l^n — 0) 
dXk+lSZi.n ■ ■ ■ SZk^n 

(68) 

If Xk+i{Qk,n, Zk+i,n = 0) = 0, take 01 into the dSD: 

_ Tf r P{Xk+l — {Xk+l} ,Qk.n, Zk+l^n ~ 9) 

l.k+i,n-e,e, J ...J dx^^^sZi,^---SZk,n 

= e^e^f (Pr [Qk.n, Zk+l,n = 0) - Pk+l,n) 

T3* 

(69) 

If i'fe+i(efe,„, Zk+i,n = 0) = Xfc+i, then the (|68ll reduces 



to: 



^k+i.n - epej/ • • •/ 



p{Xk+l = 0, 6fe,„, Zk+l^n = ^) 
SZi.n ■ ■ ■ SZk.n 



(Xfc+l - Xfc+l) (xfe+1 - Xfc+l) p{Xk+l = {Xfc+l} 



0/c,n, Zk+l.n— 9)d:X.k+lSZi.n ' ' ' SZk^n 

> e^e^ pk+i.n + 

Qk''Qk'E{Fk} [3k + E{FlQ^'Fk}r'E{Fk}Q^' 

Pr(efe,„,Zfe+i = 0) 

= epeg pfe+i_„ + [Jfc+i.n] * Pr {Qk,n,Zk+i = 0) 

~ '^k+l,n 

(70) 

The lower bound will be the minimum of the bounds of 
and (|70ll (the proof detailed to one step MSE in ifTSl ). 



B. Proof of Proposition 2 

Erom the definition of conditional probability, the relation- 
ship between Pr(0fc+i.„) and Pr(0fc.„) is as follow: 

Pr {Qk+l,n) = Pr (efc,„) p {Zk+l.n \Qk,n ) (71) 
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It is obvious that we should determine the recursion of 

P {Zk+l.n |0/£,n )■ 

When the sequence number n is in the range that 2*^+2'''^^ + 
1 < n < 2*^+^ the measurement Zk,n ^ Zk+i.n ^■ 
It means that the target exists at time step k, still exists at 
time step k + 1, and is observed. From the dynamical and 
measurement model, we can see that 



p {Zk+l,n ^ \ek,n-2^ ) = rPd, 2^= + 2'=-! 

and 



1 < n < 2'=+! 

(72) 



piZ. 



k+l,n 



Ok^n ) = 1 - rPd, 2''-^ + 1 < n < 2*= (73) 



When the measurement Zfe „ ^ 0, the state of the target 
is uncertain. Furthermore, it is hard to determine the state at 
time step k + 1 and whether the measurement is empty or not. 
Because there is the relationship that: p {Zk+i.n ^ \&k,n ) ~ 
1 — p{Zk+i^n = 0|0fe.n)- Therefore, the key is to get the 
recursive form of p{Zk+i,n = \&k,n)- This is discussed in 
the proposition 4. 

C. Proof of Proposition 3 

As in lfT4l . the recursive Bayes filter for RFS tracking 
system: 

piXk\&k-l,n) ^ J fiXk\Xk-l)piXk-l\ek-l^n)6Xk-l 

(74) 

/V in 7 \ 9{Zk,n \Xk)piXk |6fc-i,K) 
p (Afe \i^k-l,n, Zk,n ) = r;^ r (75) 

P \^k,n |C)fc-l,n ) 
p{Zk,n |6fc-l,n) = j g{Zk.n \Xk) p {Xk \Qk-ljn) 5Xk 

(76) 

Extending (|76] | by the measurement model ( l28T l and 



p{Zk.7i = |6fe-l,ri) 

= g{Zk,n='i\Xk=9)p{Xk=<J)\ek-l,n) 

+ J g{Zk,n = = {xk})p{Xk = {xfe} |efe_i,„)£ixfc 

= p (Xfe = |efe_i.„ ) + (1 - p^) (1 - p (Xk = |efc_i,„ )) 

= il-Pd)+PdPiXk=9\ek-l,n) 

(77) 

Hence dTST l reduces to: 

==0|efc_i.„,Zfc,„ = 0) 

_ g(Zfc,„=0|Xfc^0 )p(Xfc^0|efc_i,„ ) 

- p(Zfe,„=0|efe-i.„) 
_ p(Xfc^0|efc_i,„) (78) 
~ (i-Pd)+PdP(Xfc=0|efc_i,„) 
^ p(Zfc.„=0|efc_i,„)-(i-Pd) 
PdP(2fc,,x=0|ek_i,„) 

Therefore, the recursion of „ is given following: 

Pk+l,n 

= I ' " I P{Xk+l = 0, Qk,n, Zk+l,n = %)5Zi^n ' ' • SZk,n 

— f...f P{Xk+l ~ 9\Qk,n, Zk+l,n = 9) * 

p{Qk,n, Zk+l.n ^ 9) SZi^n ■■■ SZk.n 
p(Zfc + l.„^0|e;..„)-(l-Prf) /;7 -(/IIP) A 

= /.../ p,p(z,+i,„=0|e..„ ) P {Zk+i,n - W |Ui,„ j 
Pr(Ofc,„)(5Zi,„---(5Zfc,„ 

= Pr(e.,„) / • . . / £(£i.±i.*.U^(i^5Zi,„ . . . SZk.n 

(79) 



D. Proof of Proposition 4 

According to the dynamical model (IZST i and (l26T l. ( f74] l turns 
into: 

p(Xfe = 0|efc_i,„) 

= p(Xfc -0|Xfc_i -0)p(Xfe_i = 0|efe_i,„) + 

=0|Xfe_i = {xfe_i})* 
p{Xk~i = {xfe_i} |efc-i.„)c?xfc_i (80) 
= rp(Xfe_i =0|efe-i,„) + 
(l-r)(l-p(Xfe_i -0|efe-i.„)) 
- (1 - r) + (2r - l)p = \ek-i,n ) 

According to (|77]i, ^ and ([HOll, if 1 < n < 2'=-\ the 
conditional probability of Zk+i,n = is derived as follow: 

p{Zk+l,n = |6fe,n) 

= P (Zk+l.n = |Ba:_i,„, Zk.n = 0) 

= (l-Pd)+Prfp(Xfe+i =0|efc_i,„,Zfc.„ =0) 

= 1 - rPrf + (2r - l)p{Xk = |efc_i,„, Zk.n = 0) 



(81) 



VIII. Conclusion 



In this paper, a performance bound for dynamic estimation 
and filtering problem, in the framework of finite set statistics, 
is presented for the first time. This bound is recursion, and 
hence it is significant for performance evaluation of tracking 
systems. In addition, the case of P^ < 1 is taken into account, 
which makes this bound realistically. Moreover, this bound 
shows the influence of the uncertainty of target existence. 

The discussion and numerical examples show that our bound 
can obtain the enumeration PCRLB, which is the true bound 
in the case of P^ < 1 and PpA = 0, when the target remains 
form the beginning to the end. Furthermore, for some targets 
of high uncertainty, which may appear or disappearance with 
certain probability, our bound is more accurate and reasonable 
than enumeration PCRLB. To this situation, by considering 
whether the state set of the target is empty or not, the bound 
calculated in this paper is more general than previous bounds 
in the framework of random vector statistics. 
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